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Quantized, compact graphs were shown to be excellent paradigms for quantum chaos in bounded 
systems. Connecting them with leads to infinity we show that they display all the features which 
characterize scattering systems with an underlying classical chaotic dynamics. We derive exact ex- 
pressions for the scattering matrix, and an exact trace formula for the density of resonances, in terms 
of classical orbits, analogous to the semiclassical theory of chaotic scattering. A statistical analysis 
of the cross sections and resonance parameters compares well with the predictions of Random Ma- 
trix Theory. Hence, this system is proposed as a convenient tool to study the generic behavior of 
chaotic scattering systems, and their semiclassical description. 
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Quantum graphs provide a very useful tool to study 
bounded quantum systems which are chaotic in the clas- 
sical limit [|l|. Here, by attaching infinite leads we turn 
the compact graphs into scattering systems, and show 
that they display chaotic scattering PJq], a phenomenon 
with applications in many fields, ranging from nuclear 
Q, atomic [g and molecular physics, to mesoscop- 
ics [|| and classical wave scattering 0. We express the 
quantum scattering matrix and the trace formula for the 
density of its resonances in terms of the orbits of the un- 
derlying classical scattering system. These expressions 
are the exact analogues of the corresponding semiclas- 
sical approximations available in the theory of chaotic 
scattering 0Jq;0| . With these tools we analyze the distri- 
bution of resonances and the statistics of the fluctuating 
scattering amplitudes and cross sections. We show that 
graphs provide new insight on the connection between 
Random Matrix Theory (RMT) and chaotic scattering. 
Moreover, we illustrate their advantages as versatile and 
convenient tools for numerical studies of chaotic scatter- 
ing. 

Consider first a bounded graph Q. It consists of V 
vertices connected by bonds. The valency Vi of a vertex 
i is the number of bonds which emanate from it, and 
we allow only a single bond between any two vertices. 
The total number of bonds is i? = ^ J2i=i '^i- ^^'^ bond 
connecting the vertices i and j is denoted by b = {i,j). 
We shall also distinguish between the directions on the 
bond using b = (j, i) to denote the reverse direction. The 
length of the bonds are denoted by if, and we shall hence- 
forth assume that they are rationally independent. In the 
directed-bond notation Lb — Lf^. The scattering graph Q 
is obtained by adding infinite leads at the vertices of Q, 
changing their valency to Vi = Vi + 1. The leads are dis- 
tinguished by the index i of the vertex to which they are 
attached. 

The Schrodinger operator on the graph consists of the 
one dimensional Laplacian {—idx — Ab) supplemented 



by boundary conditions on the vertices Q. The "vector 
potentials" Ai, = —Af^ are introduced to break time re- 
versal symmetry, and their value may vary from bond to 
bond. On each of the bonds b or leads i the wave function 
is expressed in terms of counter propagating waves with 
a wave number k: 

On the bonds : i/'h = a6e^('^+^^)^'' + a.-e'^'^+^o^^-^"-"^'') 



On the leads : i/'i — li^ 
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where the coordinate Xb on the bond b = {i,j) takes the 
value (Lb) at the vertex i (j) while Xi measures the 
distance from the vertex along the lead i. 

The amplitudes Ob^Oj^ on the bonds and Ii,Oi on the 
leads are determined by matching conditions at the ver- 
tices. They are expressed in terms of the Vi x Vi ver- 
tex scattering matrices S^- ', , where j, j go over all the 
Vi bonds and the lead which emanate from i. The E^'^ 
are symmetric unitary matrices, which guarantee current 
conservation at each vertex by requiring 
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where Cj,,i = aj^^ie'^'^'^^'•^'il'>^^'•'■•'lK Above, the vertex 
scattering matrix S'*^ was written explicitly in terms of 
the vertex refiection amplitude p^^ , the lead-bond trans- 
mission amplitudes {tj }, and the Vi x Vi bond-bond 

transition matrix cr;'j,, which is sub unitary (| deta^'-*! < 
1), due to the coupling to the leads. As an example, 
for u-regnlar graphs {vi = wVi) with Neumann matching 
conditions on the vertices, ||l|]: 
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Combining the equations (g) for all the vertices, we ob- 
tain the V X V scattering matrix S'^^^ which relates the 
outgoing and incoming amplitudes on the leads, 



s^V^k^P^'^ 
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Here, Df^gj){k;A) = exp(i(fc + A(s j-))L(^j)) is a diago- 
nal matrix in the 2B space of directed bonds. The ma- 
trix S{k; A) = D{k\ A)R propagates the wave functions: 
R assigns a scattering amplitude for transitions between 
connected directed bonds: R(i^r),{s,j) — 5r,s^ij] D{k;A) 
provides the phase due to free propagation. The matrix R 
is sub-unitary, since |deti?| — Y[i=i Ideto-^'^j < 1. The 
scattering matrix (0) is interpreted in the following way. 
The prompt reflection at the entrance vertex induces a 
"direct" component. The "chaotic" component starts by 
a transmission from the incoming lead i to the bonds 
(i,r) with transmission amplitudes r^ . Multiple scat- 
tering is induced by (/ - S{k;A))-'^ = J2n=oS''{k;A): 
The wave gains a phase Q>'ik+Ab)Lb j^j. g^ch bond it tra- 

verses, and a scattering amplitude ct^.s at each vertex, 
until it is transmitted from the bond (s, j) to the lead j 
with an amplitude r^ . Explicitly, 
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where Vi^j is the set of the trajectories on Q which lead 
from I to j. Bp is the amplitude corresponding to a path 
p whose length and directed length are Ip = J2bep ^b and 
bp = J2bev ^bAb respectively. The scattering amplitude 
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Si is a sum of a large number of partial amplitudes, 
whose complex interference brings about the typical ir- 
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as a function of k. 



regular fluctuations of jS- ^ 

The resonances are the (complex) zeros of 



Zg{k)^det(l^S{k;A)') 
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The eigenvalues of i^ are in the unit circle, and therefore 
the resonances appear in the lower half of the complex 
k plane. Denoting the eigenvalues of S'^^^(fc) by e'^'^^'"'\ 
det S^^^k) = exp[ie(fc)] = exp[iJ2r^0r{k)] is derived 
from (0) by standard manipulations pO[ , giving 

e(fc) - 6(0) = -22:mlogdet(/ - S{k; A)) + Ck . (7) 

C — 2J2b=i^b is twice the total length of the bonds 
of Q. The resonance density dji{k) (which is propor- 
tional to the Wigner delay time) |^,0 is given by 
dii{k) = -^ ~^f^~' . It assigns to each resonance a nor- 
malized Lorentzian centered at 5Re(A:„) with a width 2 
'^m{kn). Hence, 



dnik) = — £+ -7^e V V nplprA;e'^'-''+''''> 



(8) 



n=lpG-p„ 



where the sum is over the set Vn of primitive periodic 
orbits whose period Up is a divisor of n, with r = n/up. 
Ip and bp are the length and the directed length, respec- 
tively, and the amplitudes Ap are the products of the 
bond-bond scattering amplitudes a\ ^, along the prim- 
itive loops. The mean resonance spacing is given by 
A = 27r/£. (H) is an exact trace formula for the reso- 
nance density. 
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FIG. 1. (a) The 5000 resonances of a single realization of a pen- 
tagon with A 7^ 0. The solid line marks the position of the gap 
Tgap. (b) The distribution •p(7). The solid line is the RMT pre- 
diction |13| . The difference 'P{-i) — VcuEil) is shown in the inset. 

The classical dynamics associated with Q can be easily 
defined on the bonds, but not on the vertices which are 
singular points. However, a Liouville description is con- 
structed (see 0] ) by considering the evolution of a phase- 
space density over the 2B dimensional space of directed 
bonds. The corresponding evolution operator consists 
of the transition probabilities Uh.b' between connected 
bonds 6,6', taken from the corresponding quantum evo- 
lution operator, Ub,b' — \Rb.b'\'^ ■ Due to scattering to the 
leads '^y Ubb' < 1, and the phase-space measure is not 
preserved, but rather, decays in time. Let pb{n) denote 
the probability to occupy the bond 6 at the (topological) 
time n. The probability to remain on Q is 
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Ubb'Pb'{n-l)^e-''^''PiO) (9) 



where exp(— Fc/) is the largest eigenvalue of the "leaky" 
evolution operator Ubb'- For the w-regular graph (0), the 
probability to leak to the lead per time step is r^, hence, 
Tci ~ (2/(1 -I- u))^. The set of trapped trajectories whose 
occupancy decays exponentially in time is the analogue 



of the strange repeller in generic Hamiltonian systems 
displaying "chaotic scattering" . 

The formaHsm above can be easily modified for graphs 
where not all the vertices are attached to leads. If Z, is not 

in the definition 
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attached, one has to set p'-' — -, .^ 

of S^'^. The dimension of 5(y) is changed accordingly. 

For generic graphs, the eigenvalues of the S matrix 
are strictly inside the unit circle so the resonance widths 
r„ = —2Xm{kn), are excluded from the domain r„ > 

Tgap = -21og(|A„iaa:|)/imaa:, whcrC Xmax is the largest 

eigenvalue of S{0; A) and L^ax is the longest bond. The 
existence of a gap - typical for chaotic scattering p] from 
sufficiently open scatterers- is apparent in Fig. la. The 
widths are scaled by the mean spacing A between res- 
onances i.e. 7„ = -^ so that (7)^ determines whether 
the resonances are overlapping ((7)^ > 1) or isolated 
{{l)k < 1) ( {■)k denotes spectral averaging). The dis- 
tribution of {7„}'s is shown in Fig. lb together with the 
predictions of RMT for the CUE ensemble ||l^ . In spite 
of the general good agreement, it deviates systemati- 
cally from the numerical result (see inset), and it does 
not reproduce the sharp gap in the width spectrum. 
The cross section fluctuations depend crucially on (7)^, 
which can be approximated by the classical decay con- 
stant (7)^; = 7ci (see Fig. 3b). For the v regular graphs 
discussed above, 7c/ ~ ^ j^ j^ . Changing v and V we 
can control the degree of overlap allowing to test various 
phenomena as will be described below. 



Knit) ^Jdx e'^^^'-'iMk + |) Mk - |)), 
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where dii{k) is the oscillatory part of dnik). Substituting 
(^), we flnd that the value of KR{t) equals the squared 
sum of amplitudes Ap of the periodic orbits of length 
rip — tC. A similar sum contributes to the spectral from 
factor of the compact graph Q |l|. The corresponding 
amplitudes are different due to the fact that Ap includes 
also the information about the escape of flux to the leads. 
Assuming that all periodic orbits decay at the same rate, 
one would expect K^it) « i^(i)e "'''='', where K{t) is the 
form factor for the compact system fl^ . This simple ap- 
proximation is checked in the inset of Fig. 2 (see dashed 
line) and it is shown to reproduce the numerical data 
rather well in the domain t < 5. The asymptotic decay 
is dominated by the resonances which are nearest to the 
gap, and it cannot be captured by the crude argument 
presented above. For generic graph, Kn decays expo- 
nentially but with a rate given by ^as = 7gap (the best 
flt, indicated in Fig. 2 by the dashed line, give 'jas that 
deviates by 30% from jgap)- For the graph with Neu- 
mann boundary conditions, jgap — and one expects an 
asymptotic power-law decay. (The corresponding dashed 
line in Fig. 2 shows i~^). 

Another signature of overlapping resonances are the 
Ericson fluctuations observed in the k dependence of the 
scattering cross-sections. They are one of the prominent 
features which characterize generic chaotic scattering, in 
the semiclassical limit. A convenient measure for Ericson 
fluctuations is the autocorrelation function 
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FIG. 2. The form factor Knit) for a pentagon with generic 
and Neumann boundary conditions with the same mean resonance 
width (7) J. and A 7^ 0. The data were averaged over 5000 spectral 
intervals and smoothed on small t intervals. In the inset we show 
Kii{t) for small times. Solid line correspond to the numerical data 
for the pentagon with generic boundary conditions while dashed 
line is the approximant Kjiit) fs; K{t)exp{—^c^it). 

We shall use (H) to study the resonance correlation 
function in terms of its form-factor 



C(x;^) = 
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lik-p)k (11) 



where Aj = j„ 



1. Substituting (|) in (H] 



we split the sum over trajectories into two distinct parts: 
the contributions of short trajectories are computed ex- 
plicitly by following the multiple scattering expansion up 
to trajectories of length Imax- The contribution of longer 
orbits are approximated by using the diagonal approxi- 
mation, which results in a Lorentzian with a width jEr, 
expected to be well approximated by "fci- Including ex- 
plicitly up to 71 = 3 scatterings we get, 
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where the constant G is determined by the normalization 
condition C(x = 0; i^) = 1. The interplay between the 
contributions of long and short periodic orbits is shown in 
Fig. 3a. For overlapping resonances, the autocorrelation 
function is well reproduced by the Lorentzian expected 
from the standard theory of Ericson fluctuations. The 



other case corresponds to isolated resonances where the 
contributions of short paths are clearly seen. From each 
of the various statistical measures of the resonance den- 
sity and the cross sections fluctuations discussed above, 
we extracted the effective average 7, which would fit best 
the numerical data. In Fig. 3b we compare these numer- 
ical values, with the classical expectation, and the pre- 
dictions of RMT Q. The results justify the use of the 
classical estimate for the computation of these quantities 
especially in the limit y — > 00 for fixed v jY (which is 
the analogue of the semiclassical limit). In this limit, the 
RMT and and the classical estimate coincide. 
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To investigate further the statistical properties of the 
S^' matrix, we study the distribution of scaled partial 
Wigner delay times T = ^ ai ■ The resulting distri- 
bution for various graphs with A — '^ and A j^ are 
shown in Figs. 4a, b respectively, together with the pre- 
dictions of RMT |T^. An overall agreement is evident. 
Deviations appear at the short time regime (i.e. short 
orbits), during which the "chaotic" component due to 
multiple scattering is not yet fully developed Q . 

To summarize, we presented analytical and numeri- 
cal results, on the basis of which, we propose quantum 
graphs as a model for the study of quantum chaotic scat- 
tering. Their simplicity enable us to get new understand- 
ing on the subject. Because of lack of space we defer the 
discussion of other results and further comparisons with 
RMT to a later pubhcation [|5|. 
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FIG. 3. Regular graphs with "Neumann" boundary conditions: 

(a) The real and the imaginary part of <7(x, ;^ = 1) for isolated 
(o) and overlapping (O) resonances. The solid lines correspond to 
the theoretical expression (ttS); (b) Mean resonance width (7)^, 
autocorrelation width 7£;ri the classical expectation 7^; and the 
RMT prediction ]13[, vs. V for constant valency v = 14. 
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FIG. 4. The distribution of the scaled partial delay times T 
for various graphs with "Neumann" boundary conditions. The 
dashed lines correspond to the RMT expectation ||l3|. (a) One 
channel and yl = 0; (b) V channels and A ^ 0. 
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